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The (2+l)D Georgi-Glashow (or Polyakov) model with the additional fundamental massless quarks 
is explored at finite temperature. In the case of vanishing Yukawa coupling, it is demonstrated 
that the interaction of a monopole and an antimonopole in the molecule via quark zero modes 
leads to the decrease of the Berezinsky-Kosterlitz-Thouless critical temperature when the number 
of quark flavors is equal to one. If the number of flavors becomes larger, monopoles are shown to 
exist only in the molecular phase at any temperatures exceeding a certain exponentially small one. 
This means that for such a number of flavors and at such temperatures, no fundamental matter 
f-*) ■ can be confined by means of the monopole mechanism. 

3D Georgi-Glashow model (else called the Polyakov model) is known to be one of the eldest 
and the most famous examples of theories allowing for an analytical description of confinement [jjj . 
However, the phase structure of this model at finite temperature has been addressed only recently. 
Namely, first in Ref. it has been shown that at the temperature T c = g 2 /2n the weakly coupled 
monopole plasma in this model undergoes the Berezinsky-Kosterlitz-Thouless (BKT) || phase 
transition into the molecular phase. Then, in Ref. Q], it has been shown that approximately at the 
twice smaller temperature, there occurs another phase transition associated to the deconfinement 
of W-bosons. 

In this paper, we shall be interested in the finite-temperature properties of the monopole en- 
semble, rather than the ensemble of W-bosons. Because of that, let us first discuss in some more 
details the nature of the above-mentioned BKT phase transition. At high enough temperature, 
one can apply the idea of dimensional reduction. The dimensionally-reduced theory is then the 2D 
XY-model, but with the temperature-depending strength of the monopole-antimonopole (MM) 
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interaction. Due to this fact, the phase structure of the model becomes reversed with respect 
to that of the usual 2D XY-model. Namely, at the temperatures below T c , monopoles exist in 
the plasma phase, that leads to the confinement of fundamental matter [Q, [|. At T > T c , the 
vacuum state is the molecular gas of bound MM— pairs, and consequently fundamental quarks are 
deconfined [|| . The analogy with the 2D XY-model established in Ref. El is that spin waves of the 
2D XY-model correspond to the free photons of the Polyakov model, while vortices correspond to 
magnetic monopoles. 

Let us briefly discuss the BKT phase transition, occuring at T = T c , in the language of the 
2D XY-model. At T < T c , the spectrum of the model is dominated by massless spin waves, and 
the periodicity of the angular variable is unimportant in this phase. The spin waves are unable 
to disorder the spin-spin correlation functions, and those decrease at large distances by a certain 
power law. On the contrary, at T > T c , the periodicity of the angular variable becomes important. 
This leads to the appearance of topological singularities (vortices) of the angular variable, which, 
contrary to spin waves, have nonvanishing winding numbers. Such vortices condense and disorder 
the spin-spin correlation functions, so that those start decreasing exponentially with the distance. 
Thus, the nature of the BKT phase transition is the condensation of vortices at T > T c . In another 
words, at T > T c , there exist free vortices, which mix in the ground state (vortex condensate) of 
indefinite global vorticity. Contrary to that, at T < T c , free vortices cannot exist, but they rather 
mutually couple into bound states of vortex-antivortex pairs. Such vortex-antivortex molecules are 
small-sized short-living (virtual) objects. Their dipole-type fields are short-ranged and therefore 
cannot disorder significantly the spin-spin correlation functions. However, when the temperature 
starts rising, the sizes of these molecules increase, until at T = T c they diverge, that corresponds 
to the dissociation of the molecules into pairs. Therefore, coming back to the Polyakov model, one 
of the methods (which will be employed below) to determine there the critical temperature of the 
BKT phase transition is to evaluate the mean squared separation in the MM— molecule and find 
the temperature at which it starts diverging. 

In this paper, we shall consider the extension of the Polyakov model by the fundamental dy- 
namical quarks, which are supposed to be massless. As it will be demonstrated, quark zero modes 
in the monopole field lead to the additional attraction between a monopole and an antimonopole 
in the molecule at high temperatures. In particular, when the number of these modes (equal to 
the number of massless flavors) is sufficiently large, the molecule shrinks so that its size becomes 
of the order of the inverse W-boson mass. Another factor which governs the size of the molecule is 
the characteristic range of localization of zero modes. Namely, it can be shown that the stronger 
zero modes are localized in the vicinity of the monopole center, the smaller molecular size is. In 
this paper, we shall consider the case when the Yukawa coupling vanishes, and originally massless 
quarks do not acquire any mass. This means that zero modes are maximally delocalized. Such 
a weakness of the quark-mediated interaction of monopoles opens a possibility for molecules to 
undergo eventually the phase transition into the plasma phase. However, this will be shown to 
occur only provided that the number of flavors is equal to one, whereas at any larger number 
of flavors, the respective critical temperature becomes exponentially small. This means that the 
interaction mediated by such a number of zero modes is already strong enough to maintain the 
molecular phase at any temperature larger than that one. 

Let us start our analysis with considering the Euclidean action of the Polyakov model extended 
by the fundamental, originally massless quarks. [Note that such a model can be viewed as the 
(2+l)-QCD with the additional adjoint Higgs field.] We shall first consider the general case with 
the nonvanishing Yukawa coupling, by means of which quarks acquire a certain mass. The action 
under discussion then reads 



s = Jd 3 x \ (i^) 2 + \ (D M #«) a + \ (V) 2 - v 2 ) 2 - (jD + h^^j j> 



(1) 
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Here. 



and V> = ipi (3 with the Euclidean Dirac matrices defined as 7 = —ifia, where 
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Next, in 3D, the electric coupling <?, the Yukawa coupling h, and the vacuum expectation value of 
the Higgs field r\ have the dimensionality [mass] 1 / 2 . The Higgs coupling A has the dimensionality 
[mass]. The masses of the W- and Higgs bosons are large compared to g 2 in the standard pertur- 
bative (else called weak-coupling) regime g <C r\ and read: mw = OT, tuh = rjV^X. The inequality 
g <C r\ is necessary to ensure the spontaneous symmetry breaking from SU(2) to U(l). Note also 
that for the sake of simplicity, we omit the summation over the flavor indices, but consider the 
general case with an arbitrary number of flavors. 

One can further see that the Dirac equation in the field of the third isotopic component of the 
't Hooft-Polyakov monopolc decomposes into two equations for the components of the SU(2)- 
doublet tp. The masses of these components stemming from such equations are equal to each other 
and read m q = hrj/2. Next, the Dirac equation in the full monopole potential has been shown (8| 
to possess the zero mode, whose asymptotic behavior at r = \x | 3> Tn7 has the following form: 



Xtn = { s t s n - s v s+) , Xun = °- ( 2 ) 



Here, Xn are the upper and the lower components of the mode, i.e. ip = (*-), next n = 1,2 is 

the isotopic index, v = 1,2 is the Dirac index, s + = (J), s~ = (^j, and Af is the normalization 
constant. 

It is a well known fact that in 3D, the 't Hooft-Polyakov monopole is actually an instanton Jl], |) . 
Therefore, we can use the results of Ref. Q on the quark contribution to the effective action of 
the instanton-antiinstanton molecule in QCD. Let us thus recapitulate the analysis of Ref. |9) 
adapting it to our model. To this end, we fix the gauge <& a = 7]S a3 and define the analogue of the 
free propagator So by the relation Sq 1 = —i (^jd + m q T 3 ^j. Next, we define the propagator Sm 

in the field of a monopole located at the origin, A aM [Af M — > e a% ^x^ / (gr 2 ) at r ^ m^], by the 
formula S^/ = Sq — g^ 1 ^A a ' M . Obviously, the propagator in the field of an antimonopole 
located at a certain point R, A aAI (x) = —A aAI (^x — Rj, is defined by the equation for S^/ 

with the replacement A aM — > A aM . Finally, one can consider the molecule made out of these 
monopole and antimonopole and define the total propagator S in the field of such a molecule, 
A a = A aM + A aM , by means of the equation for S M X with A aM replaced by A a . 

One can further introduce the notation \ipn), n — 0,1,2,..., for the eigenfunctions of the 
operator —tfD defined at the field of the molecule, namely —tfD \tp n ) = X n \ip n ), where Ao = 0. 
This yields the following formal spectral representation for the total propagator S: 

S(x,y) = f ^f )HM y )l . 

^ A„ - lITlqT 6 
n— y 

Next, it is convenient to employ the mean-field approximation, according to which zero modes 
dominate in the quark propagator, i.e., 
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S (x, y) 



|^,(f)) (Mv)\ 



+ S (x, y) . 



(3) 



Indeed, this approximation is valid, since in the weak-coupling regime the monopole sizes, equal 
to 7fl>yy , are much smaller than the average distance in the MM— plasma. This average distance 
has an order of magnitude £ -1 / 3 (see e.g. Ref. jl0) for a discussion). Here, C ^ e~ 47rmwe / s2 
stands for the so-called monopole fugacity, which has the dimensionality [mass] 3 , and e ~ 1 is 
a certain dimensionless function of (rnu/mw)- Obviously, £ is exponentially small in the weak- 
coupling regime under study. The approximation (^|) remains valid for the molecular phase near 
the phase transition (i.e. when the temperature approaches the critical one from above), we shall 
be interested in. That is merely because in this regime, molecules become very much inflated being 
about to dissociate. 

Within the notations adapted, one now has S = (S M + — Sq 1 ) = SfjS^ 1 Sm, where 



S — So — (Sm — So) S 1 (S M — So) — So 



-im q T° 



and 



M 



tpQ 1 ) are the zer0 modes of the operator —ijD defined at the field of a monopole and 



an antimonopole, respectively. Denoting further a = i^ip^ 
to see by the definition of the zero mode that a = ( ip, 



aM 



ipo /, it is straightforward 



'-'() 



—i'jd 



^) = (^\So 1 \^)- 



This 



yields S = S + (a*/m 2 ) 



where the star stands for the complex conjugation, and 



therefore detS = 1 + (\a\/m q 

ln[detS'- 1 /det% 1 
r = const + Nf In [rd z q 



detSo- Finally, defining the desired effective action as T 



we obtain for it in the general case with Nf flavors the following expression: 
| 2 ). The constant in this formula, standing for the sum of effective 



actions defined at the monopole and at the antimonopole, cancels out in the normalized expression 
for the mean squared separation in the MM— molecule. 

Let us further set h equal to zero, and so m q is equal to zero as well. Notice first of all that 
although in this case the direct Yukawa interaction of the Higgs bosons with quarks is absent, they 
keep interacting with each other via the gauge field. Owing to this fact, the problem of finding 
a quark zero mode in the monopole field is still valid 0. The dependence of the absolute value of 
the matrix element a on the distance R between a monopole and an antimonopole can now be 



J d 3 r/ 



R 



-4ir ln(/ii?), where /x stands 



straightforwardly found. Indeed, we have \a\ cx 
for the IR cutoff. 

Now we switch on the temperature T = /S^ 1 , so that all the bosonic (fermionic) fields should be 
supplied with the periodic (antiperiodic) boundary conditions in the temporal direction, with the 
period equal to (3. The magnetic-field lines of a single monopole thus cannot cross the boundary 
of the one-period region and should go parallel to this boundary at the distances larger than f3. 
Therefore, monopoles separated by such distances interact via the 2D Coulomb law, rather than 

1 Note that according to Eq. (^|) this mode will be non-normalizable in the sense of a discrete spectrum. However, 
in the gapless case m q = under discussion, the zero mode, which lies exactly on the border of the two contiguous 
Dirac seas, should clearly be treated not as an isolated state of a discrete spectrum, but rather as a state of a 
continuum spectrum. (A similar treatment of the zero mode of a massless left-handed neutrino on electroweak Z- 
strings has been discussed in Ref. Wm.) This means that it should be understood as follows: \ipo(x)) ~ lim (e lpr /r), 

where p = \p\. Once being considered in this way, zero modes are normalizable by the standard condition of 

oo 

normalization of the radial parts of spherical waves, R p i, which reads Jl^ ] J drr 2 R p nR p i = 2ir5(p' — p). 
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the 3D one. Recalling that the average distance between monopoles is of the order of C -1 ^ 3 ) we 
conclude that at T > C 1 / 3 , the monopole ensemble becomes two-dimensional (see e.g. Ref. (2| 
for a detailed discussion of the dimensional reduction in the Polyakov model). However, at the 
temperatures below the exponentially small one, C 1 / 3 , monopoles keep interacting by the usual 
3D Coulomb law, and the monopole confinement mechanism for the fundamental matter works at 
such temperatures under any circumstances. 

We are now in the position to explore a possible modification of the standard BKT critical 
temperature S T c = g 2 /2n due to the zero-mode mediated interaction. As it was discussed above, 
this can be done upon the evaluation of the mean squared separation in the MM— molecule and 
further finding the temperature below which it starts diverging. In this way we should take into 
account that in the dimensionally-reduced theory, the usual Coulomb interaction of monopoles ^| 

+ 00 —1/2 

R^ 1 = (7^ 2 + {ft 12 ) 2 ) goes over into — 2T\n(fj,7l), where 1Z denotes the absolute value of 



the 2D vector 1Z. This statement can be checked e.g. by virtue of the Euler - Mac Laurin formula. 
As far as the novel logarithmic interaction, ln(^i?) = ^ In \x (1Z 2 + (/3n) 2 ) , is concerned, 

71 — — OO L J 

it transforms into 



irTK + In [1 - exp(-27rTft)] - In 2. (4) 
Let us prove this statement. To this end, we employ the following formula |13|: 



00 

V 1 = — 

^ n 2 + x 2 2x 

n—l 



7r coth(7ra;) 

a; 



This yields 



X 



s^—\ 1 1 xa 1 a ^fax\ 

^ x? + O.irn. I nS 2 ~~ X + 2^ ^ n 2 4- (xal2ir\ 2 ~ 2 C ° \~2 ') 



n=~oo X 



+ {2-Knjaf x 2tt 2 ^ n 2 + ( a;a /27r) 
On the other hand, the L.H.S. of this expression can be written as 

n— — 00 \ / 

Integrating over x with the constant of integration set to zero, we get 




= az + 21n(l- e~ ax ) -21n2. 



Setting 2n/a = /j,/3 and x = [ilZ we arrive at Eq. (|]). 

Thus, the statistical weight of the quark-mediated interaction in the molecule at high tem- 
peratures has the form exp(— 2Nj In \a\) oc {itTIZ + In [1 — cxp(— 2nTlZ)] — ln2]~ 2Nf . Accounting 
for both (former) logarithmic and Coulomb interactions, we eventually arrive at the following 
expression for the mean squared separation (i 2 ) in the molecule as a function of T, g, and Nf. 

2 Without the loss of generality, we consider the molecule with the temporal component of R equal to zero. 
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/ dftTT ^ [ttT^ + In [1 - exp(-27rTft)] - In 2f 2JV/ 

/ t2\ _ ™w 

\ / _ "~55 777; • 
/ dKK [ttT^ + In [1 - exp(-27rT^)] - ln2] _2JV/ 

In this equation, we have put the lower limit of integration equal to the inverse mass of the W- 
boson, which acts as an UV cutoff. 

At large Tl, In 2 < irTTZ and |ln [1 - exp(-27rT^)] | ~ exp(-27rTft) < irTTZ. Consequently, 
we see that (£ 2 ) is finite at T > T c = (2 — Nf)g 2 /Ant, that reproduces the standard result § at 
TV/ = 0. For iV/ = 1, the plasma phase is still present at T < g 2 /Air, whereas for Nf > 2 the 
monopole ensemble may exist only in the molecular phase at any temperature larger than C 1 / 3 . 
Clearly, at Nf 3> max {l, AirT/g 2 }, (L 2 ) — > m^ 1 , which means that such a large number of zero 
modes shrinks the molecule to the minimal admissible size. Note finally that both the obtained 
critical temperature g 2 /Air and the standard one (in the absence of quarks), g 2 /2ir, are obviously 
much larger than C 1 ^ 3 : that fully validates the idea of dimensional reduction. 

In conclusion of this paper, we have found the critical temperature of the monopole BKT phase 
transition in the weak-coupling regime of the Polyakov model extended by the massless dynamical 
quarks, which interact with the Higgs boson only via the gauge field. It has been shown that for 
Nf = 1, this temperature becomes twice smaller than the one in the absence of quarks, whereas 
for Nf > 2 it becomes exponentially small, namely of the order of C 1 / 3 . The latter effect means 
that this number of quark zero modes, which strengthen the attraction of a monopole and an 
antimonopole in the molecule, becomes enough for the support of the molecular phase at any 
temperature exceeding that exponentially small one. Therefore, for Nf > 2, no fundamental 
matter (including dynamical quarks themselves) can be confined at such temperatures by means 
of the monopole mechanism. 
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